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ABSTRACT 


An Advanced Simulation of Digital Filters has been 
implemented on the 18M 360/67 computer utilizing Tektronix 
hardware and software. The program package Is appropriate 
for use by persons beginning thelr study of digital signal 
processing or for filter analysis. The ASDF orosgrams provide 
the user with an interactive method by which filter pole and 
zero locations can be manipulated. Grpahical output on both 
the Tektronix graphics screen and the Versatec plotter are 


provided to observe the effects of pole-zero movement. 
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1. INTRONUCTION 


Dieltal methods have recently begun to dominate’ the 
selection of transmission and filtering schemes. When the 
student first confronts the discrete domain, he encounters a 
series of unfamiliar concepts. While most of the discrete 
concepts have parallels in the continuous or s-domain, this 
dualism is not apparent inittally. 

The programs comprising the Advanced Simulation of 
Digital Filters (ASDF) provide the user with a convenient 
method for developing intultion in the discrete domain. 
This is achteved by allowing the user to generate digital 
filters by directly entering the ooles and zeros into a 
graphics terminal. The ASDF programs then perform the 
calculations required to display the characteristics of the 
user generated filter. 

The overriding objective of the ASDF program package Is 
to provide the user with an interactive method for modifying 
the filter pole/zero locations [tn the z-plane, as well as to 
compute the filter responses. Both time and frequency 
responses are computed and displayed graphically on the 


Tektronix 4912 graphics terminal. (Ref. 9] 
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tt. SELECTION OF ASDF ALGORITHMS 


There are several parameters which are particularly 
important in selecting the algorithms used in the ASDF 
programs. Since one of the primary objectives of these 
programs is to provide the user with clear explanstions of 
methods for computing digital filter charactertstics, direct 
or strafghtforward methods should be used. The Importance 
of developing a user's intuition for the discrete domain and 
digital filters suggests that graphical rather than 
numerical outputs would be most appropriate. The use of 
graphical outputs to characterize a filter's responses 
Impltes minimum Importance of precision in these results. 
Since the reason for computing the filter characteristics is 
to develop the user's ftntultion for the discrete domain, 
there are no pressing constraints on the selected algorithms 
with regard to computational efficiency, speed, or memory 
space utilfzatton. In fact, the final ASDF package spends 
the greatest portion of its execution time making the 
interface with the z-plane, and the form of the output easy 
to understand and to use. The greatest part of the memory 
space utilized for the ASDF programs [is consumed in the 


Interfacing processes rather than tn the methods’ by which 


the data fs stored, used, or computed. 


A. ALGORITHMS SELECTED 

Because of the above stated objectives, a simple 
algorithm was selected for computation of time domain 
responses of digital filters. Once the user enters the 
destred pote and zero locations, a direct-form-two 
difference equation representation of the filter Is 
generated [Ref. 5]. The unit sample sequence is used for 
the filter input, the difference equations are solved 


recursively in the discrete time domain and the unit sample 


response Is recorded as the output. The unlit step response 


is calculated by a similar method. 

It Ts known that the direct-form-two filter 
implementation is very sensitive to parameter variations, 
and as such cannot be utilized tn sttuations requiring 
extreme accuracy. The advantages associated with simplicity 
of understanding the algorithms, the ease of programming the 
general case, and the deemphasis of accuracy In outputs for 


marginally stable filters outweighs this disadvantage. 


B. DEFINING THE DATA STRUCTURE 

A detailed description of the Important variables in the 
ASDF programs is ziven [nm Apnendix 8. There is, however, an 
overiding data structure common to all portions of the ASDF 
package which is Important. 

The array POLZRO holds the locations of the poles and 


zeros of the filter under analysis. tn all cases, zeros are 
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stored in POLZRO(I,J) for values of | from one to ten. Pole 
locations are stored in POLZR29(1,J) for values of I from 
eleven to twenty. Only the top half plane pole of each 
complex pair is stored. There are five sections of 
POLZR0(!,JU), }.e. values of J from one to five. POLZRO(I,1) 
is the real part of. the rectangular representation of the 
ith root. POLZRO(1,2) is the imaginary part of the 
rectangular representation of the ith root. POLZRO(!,3) and 
POLZRO(!,%) are the magnitude and angle portions of the [th 
root location in polar form. POLZRO(1,5) Is the root tyne 
designator. Root type des!fgnators are: one for real zeros, 
two for complex zeros, three for real poles, and four for 
complex poles. 

The array !ROOTS(!) holds the numbers of roots of each 
type of root already entered into the system. JtROOTS(1) and 
{R00TS(3) hold the number of real zeros and poles, 
respectively. !ROOTS(2) and !ROOTS(4) hold the number of 
complex zero and pole pairs, respectively. 

The variable !ERROR is the error flags, which is” set 
whenever the user generates an error. There are numerous 
possible errors. These fnclude trying to delete roots which 
have not been entered, or adding roots which cause’ the 


maximum system order to be exceeded. 
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111. IMPLEMENTED SOLUTION 


Having determined the general form of the solution it 


remalned to generate the software Implementation. There 


were five major areas to be programmed: the interface with 


the time response, the frequency response, the 


the z-plane, 


plotting of output curves, and the supervisory routines. 


A. INTERFACE WITH THE Z-PLANE 


Stnee the method by which the user modifies the pole and 


zero locations in the z-plane is of no theoretical fnterest, 


will be brief. 


the discussion of this implementation 


Extensive coverage of the subroutines employed and the 


is reserved for Apnendix 


definitions of the variables used 


B. 


The ASDF command which allows the user to make pole-zero 


modifications in the z-plane is POLZ&9. This command 


initiates execution of a load module which ts structured as 


depicted in figure 3-1. Each of the boxes shown represents a 


on both ends indicate 


with arrow heads 


subroutine. Lines 


program control passes to the subroutine and then 


that 


to the calling routine «:th no intermediate 


dtrectly back 


branches. The main orogram, MAIN, initializes the required 


arrays, and tests to ensure that the user has’ generated a 


realizable causal filter. 


_ 
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Figure 3-1 Flowchart Of POLZRO 


The actual interaction with the z-plane Is accomplished 


In subroutine UNTCIR. The structure of UNTCIR [Is diagrammed 


in figure 3-2, The large double ended arrow indicates that 


there are many subroutines of the form CRRMMM, The 


subroutines LETTER and BOXUC generate the unit circle and 


command array. PLOT-10 software (Ref. 10] statements 


determine the location of the user positioned cursor. This 


cursor position fs adjusted with software tabs to one of 


twenty-six possible locations. Based on the X and Y 


location of the adjusted cursor position a specific 


subroutine ts selected for execution. The flashing 


alphanumerte cursor ts positioned in the selected command 


box using the adjusted cursor location as a reference. The 


subrout}]ne then reads two alphanumeric characters typed by 


the user. The combination "space", "e" indicates that the 


user wishes the selected subroutine to be executed. The 


combination "space", "x" indicates that the user wishes to 


terminate the pole zero manipulation phase of the program. 


All of the commands which change the status of the 


z-plane have names of the form CRRMMM where C [s either "A" 


for add root(s) or "D' for delete root(s). The RR Indicates 


y: root type: RZ - real zero, CZ - complex zero, RP - rea) 


! pole, CP - complex pole. The MMM stands for the mode of 


data entry. !AG stands for tnteractive graphics i.e. the 


roots will be adjusted with the cursor. POL and RCT stand 


for polar and rectangular, respectively. The general 


OLFLTR 


ERROR 


LETTER 


BOXUC 


UNTCIR PLOTUC 


Figure 3-2 


CRRMMM 


Structure Of Subroutine UNTCIR 
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structure of the CRRMMM subroutines {fs described In figure 
3~3. 

Al} CRRMMM subroutines follow the same general pattern. 
The POLZRO table Is checked to determine whether or not the 
issued command fs valid. Poles and zeros must be present 
before they can be deleted. There cannot be more than ten 
poles or ten zeros. The desired pole or zero location is 
next acquired. For subroutines adding roots! the entered 
location is. that of the root(s) to be added. For 
subroutines deleting roots the entered tlocation is any 
location near the root(s) to be deleted. Using subroutine 
LOCATE, the POLZRO table is searched for the root of the 
correct type which is nearest the specified locatton. This 
nearest root Is then highlighted. In all cases, except when 
adding roots in’ the tAG mode, the user [ts asked whether or 
not the correct location has been determined. Once the root 
location ts verified the root is added and plotted with 
subroutines STRRTS and PLOTRT. For subroutines deleting 
roots the POLZRO table is adjusted with subroutine UPDATE 
and the unit circle is redrawn. 

The user may make an error in carrying out’ the above 
steps, tI.e., by placing a pole outside the unit circle, by 
exceeding the maximum system order, or hy entering a zero 


outside ten units from the origin. In each of these cases an 


luereafter, poles and zeros will be refered to as roots. 
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ERROR 
RETURN 2 


No 


PRINT 
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Figure 3-3a Structure of CRRMMM Subroutines 
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ROOT 
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UPDATE 


RETURN 


Figure 3-3b Structure Of CRRMMM Subroutines Continued 
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STRRTS SET IERROR 


PLOTRT 


RETURN 


Figure 3-3c Structure Of CRRMMM Subroutines Continued 
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error flag is set. Subroutine ERROR handles the printing of 
messages and subroutine LST prints a Jisting of the current 


filter status before the program resumes. 


B. TIME RESPONSE 
The transfer function of a digital filter Is defined In 


equation 3.1. The Zis Pus ay, and b, will be defined ina 
subsequent section. The variable A is an arbitrary constant 


assigned by the user and is called the filter gain. 


M M N 
An (2-z.) Eb. 2°7 
i=l - rao 
H( 2) = ————_—_—_———_-# —————_——_—_———_ (3.1) 
N k 
T (z-p.) l- £ ay. Zz 
tel = k=l 


Inherent In generating the time domain response with a 
direct-form-two filter Implementation are the determination 
of the coefficients ay and b. In equation (3-1), and the 
computation of the unit sample and unit step sequences. 
Having determined the filter coefictents, the time domain 
response of the flliter ts determined by using the 
appropriate filter input sequence and recording the output 
of the filter for the desired number of points. This Is 
accomplished by recursively solving a difference equation 
which characterizes the filter for 1024 ttme points. The 
standard direct-form-two filter [mplementatlion is shown In 


figure 3-4. The zt terms represent unit time delays. This 
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figure corresponds to the difference equations: 


N 
y(nek) + f b.x(n-Ner) (3.2) 


N 
y(n) = fa 
ks r= 


1 k 

To compute the fiiter's time response, three operations 
are required: calculate the a, and ae coefficients, generate 
the [tnput sequence x(k), kz0, 1, . . . 1024, solves the 
difference equations and stores the filter output values. 
There are four central subroutines which perform these 
operations: ASUBK, BSUBR, RSPNSE, and COEFF. Subroutines 
ABSUK, BSUBR, and COEFF generate the filter coefficients 
starting with the oole and zero locations specified. 
Subroutine RSPNSE solves the difference equations by 
coordinating the other three routines, generating the 
appropriate Input sequence, and storing the output sequence 
values. 

The ASUBK and BSUBR subroutines are f[dentical [In 
structure and purpose. ASUBK computes the degree of the 
denominator polynomial as two times the number of complex 
pole pairs plus the number of real poles. A complex vector 
ts formed named POLES, which holds a listing of all the 
poles of the filter under analysis. The POLZRO array stores 
real poles in’ the form X+#j0. Only the top half olane pole 
of each complex pole pair fs stored, t.e., (X+jY), where X 
and Y are non-negative. The vector POLES has entries for 


all real poles and both poles of each complex pole pair. 
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Once the POLES vector has been generated the subroutine 


COEFF multiplies out the factors of the form (z-p,) to 
generate the coefficients of the denominator polynomial. 

The algorithm for multiplying out the denominator 
polynomial is essentially the same as would be done by hand. 
multiplication of a series of literal factors of the form 
(z-p;). The vector POLES Is oassed to the COEFF subroutine 
as the vector V. Since the denominator factors are of the 
form (z-p5), the V vector is changed to hold the values “PP; 
If the system is of order less than ten, the remaltning 
positions of V are filled with zeros. The degree of the 
denominator polynomial Is subsequently used to determine how 
many poles are actually located at the origin. 

The vectors used in the subroutine COEFF hold the 
coeffictents for the polynomial being eenerated. The 
literal powers of "z" are supressed. The vector H1 always 
contains the partial product of the K factors oreviously 
multiplied together. The vector H2 represents the old 
partial product times the "z" of the (z-p;) factor. Since 
the coeffictent of '"z" ts always one, H2 is a left shifted 
version of H1. The vector H%3 represents the partial product 
times the next pole location (p;). This next “Pp; value Is 
stored in V(1). The new H1 vector ts formed hy adding the 
H2 and H3 vectors. H1 now contains the product of the first 


K*1 factors. The values In the V vector are next shifted 


such that the next pole location to be multiplted is always 


in the V(1) posttion. These four steps are repeated ten 
times to menerate the coefficients of the complete 


"2" The formula for the transfer function 


polynomial in 
requires coefficlents of the denominator polynomial in terms 
of z~! so the order of the coefficlents is reversed. 

Figure 3-5 depicts several iterations of the algorithm 
used in the COEFF subroutine. Figure 3-5a shows the [Inittal 
conditions for each of the vectors defined in the above 
Oaragraph. Figure 3-S5c shows the values of the vectors 
after the first Iteration. Figure 3-5f shows the vector 
values after the second iteration. The third and subsequent 
iterations follow this pattern. 

In orecisely the same way as discribed above, the ay 
coefficients are senerated with the ASUBK subroutine. The 
zero locations are stored fn POLZRO exactly as are the 
poles. The generation of the vector ZEROS and the use of 
the subroutine COEFF parallels the 8SUBR subroutine exactly. 

The stgnificant variables in subroutine RSPNSE are: 
XSUBN, INPUT, OUTPUT, FLTRGN, A, and 8B. The vectors A and 8 
hold the already computed values of a, and o,. A(1) holds 
@ 3 AC2) holds a, . . . A(11) holds aj. The b, are stored 
similarly in the B vector. XSUBN is the current value of the 
input sequence. INPUT is the value of the node so labeled 
In figure 3-4. OQUTPUT jis similarly noted on figure 3-4, 
FLTRGN Is the value of the filter gain, and its assigned by 


the user. The vector Z holds the delayed INPUT values. 
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Figure 3-5a. Initial Conditions Of COEFF Vectors. 


Figure 3-5b. 
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Values Of COEFF Vectors S3efore The End Of The 


First Iteration. 
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Hl - - ; 2) 0 1 -p 
H2 : ; z 0 9) 1 D) 
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Figure 3-5c. Values Of COEFF Vectors startine The Second 


Iteration. 
V 
H1 : . 3 0 i) 1 “Pp 
H2 F ‘ 5 9 1 -p 0 
H3 ‘ . 9 f) -p #p.p 


Figure 3-Sd. Values Of COEFF Vectors before The End Of The 


Second Iteration. 
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H1 ‘ ; F 1 -(9,4n,) (pip) 


H2 : P : 1 -p 0 
H3 ‘ . P 9) -:—p (p, p,5) 


Figure 3-S5e. Values Of COEFF Vectors Starting The 


Third {fteration. 


V 
H1 : : 1 -(p, +05) (p, Pp) 
H2 ; 1 -(p,+p,) (p, pp) ») 
H3 : 0 “p +(p, #9, )p, =(P;P.P,) 


te Figure 3-5f. Values Of COEFF Vectors before The End Of The 


Third Iteration. 
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Z(1) ts no delay; 2(2) is delayed one time unit; 2Z(3) is 


delayed two time units, etc. 

The subroutine RSPNSE senerates the appropriate input 
sequence. For the unit sample response the input sequence 
is 1, 0, 0, 0, ... . For the unit step response the input 
sequence ts 1, 1, 1, 1, .. .. The appropriate value of 
XSUBN is determined by the type of input sequence. The 
variable INPUT is then computed to be the sum of XSUBN and 
the a's times the delayed INPUT values. The filter output 
is the sum of the b's times the delayed INPUT values. 
These output values are then multiplied by the filter gain, 


FLTRGN, and stored in the appropriate array: Yl for the 


unit sample response, YS for the unit step response. 


C. FREQUENCY RESPONSE 

For discrete systems, the frequency response is found by 
evaluating the system transfer function at points on the 
unit circle. Subroutine FROANCY eenerates both the phase 
and the magnitude of the transfer function. 

The Important vartables in the FRQNCY subroutine are: 
AA, 8B, THETA, Z, NNUM, and NODEN. The variables AA and BB 
are the double precision values of the A and B_ vectors 
described in’ the section on time response. THETA is the 
angular displacement ground the unit circle for which a 


value of the frequency response is calculated. Z is the 
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classical z-plane vartable: 
Zzexp(j THETA) (3.3) 


NNUM is the number of numerator terms. NDEN Is the number | 
of denominator terms. | 

The frequency response is computed for 1024 points | 
around the unit circle. All points computed are for THETA | 
greater than zero and less than 7. The values of 2Z for 
THETA greater than -n and less than zero are wneither | 
computed nor plotted since the magnitude of the transfer | 
function is an even function of THETA and the phase response 
is an odd function of THETA. 

The generation of the a, and b,, coefficients has already 
been discussed in the section on time domain response. The 
evaluation of H(z) is straightforward. The value of THETA is 
incremented through 1024 steps from zero to 7. For each 
value of THETA generated, the value of 2 is computed. The 
value of H(z) is then computed in two parts. The numerator 

r is the sum of NNUM terms of the form b, times z* for 
U r20, 1, 2,. . . NNUM-1 ‘ The denominator Is one minus the 
’ sum of (NDEN-1) terms of the form a, times 2* for 


| keels 258 35. «s. 468 NDEN-1°. Dividing numerator hy denominator 


' 1 wNuM Ts equal to the degree of the numerator plus one. 


2 NDEN equals the degree of the denominator plus one, 
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yletds H(z). The phase of H(z) is then determined using the 
principal values of the arctangent function; the results are 
stored In the vector YP, The magnitude of H(z) Is computed, 
multiplied by FLTRGN, the filter sain, and subsequently 


stored tn the vector YM. 


D. OUTPUT PLOTS 

The output plots generated by ASDF are generated by two 
different sof tware programs each being different 
implementations of the same plotting requirements. One set 
of plots Is generated using PLOT-10 software on the 
Tektronix 4012 terminal. Another set of plots utilizes the 
Versatec plotter and associated software. [Ref. 7] The 
PLOT-10 software is adaquately covered in reference 10. The 
Versatec software Is explained In Technical Note-034 in the 
computer center. 

The plots which are generated on both devices are 
essentially the same. There are five plots generated on the 
4012 terminal: unftt sample response, unit step response, 
phase response, magnitude of the transfer function, and 
magnitude of the transfer function in decibels. Both of the 
time domain plots are terminated when the output value has 
stabilized to within one percent of the maximum senerated 
system output. All 1024 points of the ohase and magnitude 


plots are plotted. 


The form of the Versatec output fs similar to that of 


the 4012. One additional olot is included. It shows the unit 


circle with poles and zeros of the filter being analyzed. 


Plotting Is terminated when the filter output stabilizes to 


within one percent of the maximum system output generated. 


On both devices, the values of unit sample response are 


subroutine PLTIMP. 


a vector Y! and plotted with 


stored in 


The unit step response values are stored in YS, and plotted 


is stored In 


PLTSTP. 


The 


with subroutine phase response 


The magnitude of the transfer function Is stored 


vector YP. 


in vector YM. The magnitude of the transfer function [fn 


in YMDB. The three vectors YP, YM, and 


dectbels is stored 


YMDB are all plotted with subroutine PLTT2F. The Versatec 


is senerated with subroutine UNTC®. 


plot of the untt circle 


The unft circle and command array are 


4912 plot of the 


generated with the subroutine UNTCIR. 


SUPERVISORY ROUTINES 


There are three supervisory routines used with the ASDF 


, programs. These supervisory routines nake the loading of 


ASDF transparent to the user. There its a substantial amount 


of flle mantpulation required to store and load the ASDF 


programs. A knowledge of the hardware and software 


for these purposes is not narticularly useful to 


available 


the user. For these reasons, the supervisory routines are 


Two of these routines are based on 


discussed only briefly. 


CP/CHS and IBM EXEC environment; the third is written in the 
(BM-360 Job Control Language. 

The JCL program moves the ASDF program from its 
permanent storage location on DISK02 to the CP/CMS disks 
OSX001. The two EXEC routines generate an appropriate 
CP/CMS work space, and man!Ipulate the ASDF text file into a 
usable form. 

The EXEC program ASDF191 resides on the user 191-P-disk. 
This routine automatically requests and formats ten 
cylInders of temporary disk space. The ASNF TEXT file fs 
moved from the OSX901 disk to the newly created temporary 
disk space 192-T. Once in place, The ASDF TEXT file name is 
changed to ASDF EXEC to generate the next file with the 
correct fille type. ASDF EXEC is then edited to remove the 
first set of data cards which hold the ASDF192 EXEC program. 
The 192-T-disk is released and the user is requested to type 
two commands to the system. The "login 192 P"™ command makes 
the 192 disk the Pedisk. The "EXEC ASDF192" command begins 
the execution of the second EXEC file. 

The meneral ASDF192 EXEC structure is depicted In figure 
3-6. ASDF192 EXEC is created by ASDF191 EXEC and resides on 
disk 192. tn a manner” sfmitar to ASDF191, ASDF192 splits 
the ASDF TEXT file, which contains all of the component 
parts of the ASDF program, tnto a sertes of smaller files. 
Once all files have been separated, titles are altered 


appropriately, text files are loaded [nto core, and load 
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COVER 


OPTIONS INTRO 


POLZRO ASDF192 RESPONSE 


EXEC HRDSCPY EXEC PUNCH 


JOBSCRD 


Figure 3-6 Structure Of ASDF192 EXEC 


modules are generated. ASODF192 then reverts to an infinite 


loop, reading and implementing the desfred user command. 


There are six user commands) available to the user: POLZRO, 


RESPONSE, EXEC PUNCH, EXEC HRDSCPY, OPTIONS, and INTRO. The 
applicability of these commands its briefly discussed i[n 
chapter four. This brief description is also presented by 


executing the OPTIONS command once ASDF has been loaded. 


yeah Ph 


IV. STRUCTURE OF ASDF COMMANDS 


The ASOF program generates the time and frequency 
responses of realizable digital filters of order ten or 
less. Chapter three describes the algorithms used for these 
computations. The procedure for executing the ASDF program 
generally consists of two distinct steps. The first step is 
to enter the poles and zeros of the filter to be analyzed. 
The second step is to perform the calculations required to 
generate and plot the desired filter responses. There are 
six commands) which may be executed within the ASDF 
environment providing the user with the ability to control 
the form of the output, to enter the poles and zeros of the 
filter, and to clarify most questions pertaining to the 
meaning and use of the zeneral commands. The six 
fundamental commands avallable to the user are: OPTIONS, 
INTRO, POLZRO, RESPONSE, EXEC PUNCH, and EXEC HRDSCPY. Only 
the POLZRO command has a_ subordinate command structure. All 
of these commands are issued by typing the command on the 
terminal keyboard in response to: *#*ASDF-RE. 

The OPTIONS command provides a listing of all commands 
avallable to the user. Each command is accompanied by a 
brief description of Its Intended use. 

The {INTRO command writes on the screen six pages of 


tnformation pertaining to the POLZRO command substructure. 
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INTRO Includes a general description of the data input 
modes, and a description of several input situations with 
the appropriate executable commands. 

The POLZRO command allows the user to modify or create a 
filter In the z-plane. All adjustment of noles and zeros 
must take place with this command. By issuing this command 
the user may add and delete poles and/or zeros as well as 
adjust the filter gain constant. The commands which are 
issued within the POLZRO environment are described in a 
later section of this chapter. 

The RESPONSE command displays the characteristics of the 
ftlter being created by the user. Five displays are 
oresented: unit sample response, unit step response, phase 
of the transfer function, magnitude of the transfer 
function, and the magnitude of the transfer function in 
decibels. 

The EXEC PUNCH command allows the user to punch a card 
deck of the filter pole and zero locations on the offline 
punch. This punched output can be nicked up i[n room [-140 
of the computer center. The cards nunched by this command 
may be subsequently reentered Into the system so that this 
particular fllter can be referenced or modified at a later 
date. 

When the user desires to enter a filter in this way, the 
punched cards are presented to one of the computer operators 


over the counter [In room [-149, The cards are punched in 
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the correct format for direct submission into the CP/CMS 
card reader. The operator subsequently reads this deck Into 
the user's virtual card reader. As the user signs onto a 
terminal and gains access to CP/CMS, the computer Includes 
an additional message in the standard login typeout showing 
that the user has a virtual reader file. This message is of 
the form: 
FILES:- 012 RDR, NO P&T, NO PUN 

When the user wishes to access this old filter data he must 
execute the CP/CMS command OFFLINE READ »# «, The OFFLINE 
READ command will read the punched card filter data from the 
user's virtual card reader, and generate a file on the 
user's disk space entitled FILE FTI1F901. Any Files with 
this name present before the OFFLINE READ command is issued 
will be over-written by the new fille. If the user does not 
wish to lose an older FILE FTI1F001, he should use the 
CP/CMS ALTER command to change the name of the old file that 
he wishes to protect. Once the user has successfully loaded 
the ASOF program, the ASDF stgnature f.e., *ASDF-2E, is 
oresented. The OFFLINE READ command may be issued any time 
that this signature appears. 

The EXEC HRDSCPY command ts avallable only when the 
SUBMIT function of CP/CMS is functioning. When the SUBMIT 
command Is not functional a message is sent to the user 
during the login procedure which states that SUBMIT is not 


available. the output from this command is ultimately 
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available in the computer center, and consists of two parts. 


The first part of the output Is the numerical listing of the 
values of the responses computed by the ASDF program. The 
EXEC HRDSCPY command is the only command which will altow 
the user access to these numerical values. The second part 
of the output is a series of Versatec plots, each similar to 
the screen presentations of the RESPONSE command, with one 
additional plot of the unit circle and pole/zero locations. 
When the EXEC HRDSCPY command Is executed the ASDF program 
will ask several questions, using the responses to generate 
a JOB card. The response to the "ENTER DESIRED JOB NAME" ts 
the title which appears on the final output. ASDF then 
submits the filter characteristics to the operating system 
for processing. The hard copy output will be found in the 


standard output boxes. 


A. SUMMARY OF DATA INPUT MODES 

There are three modes through which the poles and zeros 
may be entered into the z-plane: 1) interactive graphics, 
2) rectangular, and 3) polar. Each of these modes’ has 
particular attributes and limitations. 

1. Interactive Graphics (IAG) 

This mode was designed to enter noles and zeros 

rapidly, in cases where accuracy is not of prime imnortance. 


Tre interactive eraphics mode is characterized by the 


appearance of the graphics cursor (the intersection of two 


lines which trace across the terminal screen). The position 
of the cursor is controlled by the two thumb wheels located 
to the right of the interactive ASCII keyboard on the 
Tektronix 4012 terminal. The coverage of the interactive 
graphics cursor is limited to values of X and Y greater than 
-1.1 and less than #1.1. Zeros outside this area must be 
entered with elther the polar or rectangular commands. 
Careful adjustment of the cursor can usually locate poles or 
zeros within .002 units of the desired location. Poles must 
be located within the unit circle. 
2. Rectangular 

This mode was designed to enter from the terminal 
keyboard, poles and zeros which are known jin” rectangular 
coordinates. tt is slower than the interactive graphics 
mode because there Is a significant amount of screen 
rewriting required. The rectangular mode allows much 
®reater accuracy In the location of voles and zeros. Poles 
may not be entered on or outside the unit circle. Zeros nust 
be located within ten units of the z-plane origin. 

3. Polar 

This mode was designed to enter from the terminal 
keyboard, poles and zeros which are known in polar 
coordinates. This mode is also slower than the interactive 
sraphics mode. It allows poles and zeros to be entered to 
any permissible position in the z-plane. Poles must be 


entered within the unlt circle. Zeros must be entered within 
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ten units of the 2z-plane origin. Each volar command has a 
description of the constraints on the data to he entered 
automatically written on the screen, This mode is 
particularily appropriate when constructing all-pass linear 


phase filters. 


B. THE POLZRO COMMAND SUBSTRUCTURE 

In an attempt to maximize execution sneed, a minimum of 
information is written on the screen, thus reducing the 
rewrite time after execution of each commaid., When orogram 
execution begins, the user is presented with the unit 
cirele, and a listing of commands on the top left side of 
the screen. The commands are abbreviated in the following 
manner. The general command type or mode is specified 
“across the top of the array as INTERACTIVE, RECT, and POLAR 
(see figure 4-1). Each of these modes conforms to the above 
description. Down the left side of the array are the 


commands of the form C8R, where C is the change to be made, 


A - for adda root, 0 - for delete a root. R& is for the 
type of root to be acted upon by the command: RZ - real 
zero, CZ - complex zero, RP - real pole, and CP =- complex 


pole. As an example, the box in the center column aligned 
with the row titled DCZ stands for the command "Delete a 
Complex Zero - Rectangular Coordinates. There are two 
unlque commands listed as NUC, and LST. The New Unit Circle 


(NUC) command simply redraws the unit circle and command 
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Figure 4-1 Basic Display Format 
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array. This is appropriate when either too many entries 
have been made to the terminal, and the command array is no 
longer clear, or when noise in the line from the computer 
causes a deformed array or unlit circle to be drawn. The LST 
command provides the user with a listing of the current 
status of the fllter being analyzed, in a self explanatory 
format. 

In all cases when the unit circle and command array just 
described appear on the screen, an audible tone or bell is 
heard. This bell Indicates to the user that the cursor is 
to be positioned within one of the boxes in the command 
array. In order to execute a command, the cursor is 
positioned in the box associated with the desired command. 
A single keyboard character (any character but the "return" 
key is acceptable) Is then typed. The cursor then 
disappears. The system presents a carat pointing to the 
command that the computer [is about to execute; the flashing 
cursor blinks within the box associated with the command 
which is about to be executed. If the box so Indicated is 
the command that was selected, the user types a “space”, an 
"te" (for execute) and then "carriage return". In cases 
where a rectangular or polar mode command is to be executed, 
the screen is. then erased, and new information is orinted. 
In cases where an interactive eraphics command is to be 
executed, only the cursor reappears. This second appearance 


of the cursor is not accompanied hy a bell or audible tone. 
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This Indicates that the cursor is to be positioned within 
the unit circle. Once so ositioned, a single keyboard 
character is typed (not "return"), and, as appropriate, the 
unit circle ts redrawn to reflect the change in pole & zero 
locations. After this update of the z-plane, the circle and 
command array are redrawn with the cursor accomnanied by the 
audible tone, 

There are only two exceptions to this scenario. First, 
if the flashing cursor does not flash in the correct box, 
the wrong command will be executed. !f an incorrect box 
contains the flashing cursor, the user should not type the 


"ea" sequence. Instead, he should type two other 


"space", 
characters (e.g., "space", "space") and then "return." The 
cursor returns along with the audible tone, calling for a 
readjustment of the cursor, closer to the desired command 
box. The second exception occurs when the placement of all 
of the poles and zeros is completed. In this case, after 
aligning the cursor with any box, the sequence "space", 'x'! 
(for exit), "return'' should be typed. This will terminate 
the input, or modification, stage of the program and present 
the user with a serfles of options for the next stage of 


filter processing. At this point a display appears which 


explains the options available to the user. 
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SUMMARY OF POLZRO SUBCOMMANDS 

When the cursor appears along with an audible tone, the 
user is to position the cursor within the command array, and 
type a keyboard character (not “return"). When the flashing 
cursor appears within a box, he types ("space", "el, 
"return") to execute the specified command. "Snace" and any 
other character but "x" will return the user to the command 
table for cursor readjustment. To terminate the Input, or 
modification mode, the cursor Is positioned over any box, 
and once the flashing cursor enters the box, the characters 


"space", "x", “return are typed. 
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V. LOADING THE ASDF PACKAGE 


The orocedure for using the ASDF program consists of 
either six ‘ar seven steps depending on whether or not the 
ASDF package ia Seas tdent on OSX901. Usually the program 
sections of ASDF are stored on the T2314 disk pack 0SX901. 
Each Monday morning the 9SX901 disks are scratched so that 
the ASDF program must be moved from the permanent storage 
location on D!ISK92, to the temporary storage location on 
OSX001. The user may simply stgn onto a 4012 terminal, and 
execute three commands to load ASOF [nto his’ virtual 
machine. Once this loading process has been successfully 
completed, the commands oreviously discussed are at’ the 


user's disposal for filter analysis. 


A. LOGIN PROCEDURE 

“ith one minor exception, the standard login procedures 
(specified in Section 4.4.3 of the CP/CMS User's Manual) 
CRef. 4] are applicable when entering the CP/CMS system to 
use ASDF. Due to the size of the ASDF programs the user 
must request a 400k virtual machine when loring onto CP/CMS. 
The virtual machine is requested as follows: 7 
Togin XXXXTNN 400k 


where XXXX is the user's computer center number, T ts the 


type of user i.e. G for seeneral user, and P_ for private 


uy 
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user, and NN is” the two digit terminal number. The 400k 
request in the LOGIN ts absolutely mandatory. Without this 
request the user will be able to load the ASDF program, 
however, when attempting to execute the RESPONSE command, an 
error: "!NPUT COMMAND ERROR" will result. Should this error 
occur when the command is corrctly spelled the user has no 
alternative but to logoff and log hack on to the system 
requesting the correct size virtual machine. The remainder 
of the login procedure is identical to the examples in the 


CP/CMS user's manual. 


B. TESTING OSXNO01 FOR ASDF PROGRAMS 
The second step in loading the ASDF programs is to 

determine whether or not a copy of the programs exists on 
OSX001. The ASDF programs are permanently stored on the 
operating system DISKO2. The first time that the programs 
are used each week, the user must move them from DISK02 to 
the disk OSX101. To determine whether or not OSX901 has the 
ASDF package, [issue the command: 

TOCP ASDF191 EXEC Pl 
The system will respond: 

EXECUTION BEGINS... 

237 = OSX001 

ENTER DSNAME OR '/s! 
The correct user response is: 


F718. asdf191 
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At this time one of two responses will appear. The 


on whether or not 


particular response recefved Is dependent 


ASDF resides on OSX901. 


1. ASDF 


is on 9SX001 


ASDF programs are located on OSX001, the 


When the 


system response is: 


EXTENTS EXHAUSTED, EOF SIMULATED. 


4 BLOCKS 63 RECORDS READ 


This indicates that file ASDF 191 has 


message 


satisfactorially been transferred to the user's virtual 


the user is ready to continue with the 


machine and 


subsequent loading steps. 


2. ASDF is not on 0SX901 


When the ASDF programs are not on OSX901, the system 


response will be: 


DATA SET NOT FOUND (08/360 $213) 
ENTER DSNAME OR '/«! 


This response indicates that the user must move’ the ASDF 


programs to OSX0N01 before the tloading procedure may be 


resumed, The simplest method for accomplishing this 


transfer is to generate a file on CP/CMS which exactly 


duplicates the Job Control Language program whitch follows. 


To tenerate the file which transfers ASDF, logon to 


fe CP/CMS. Once the system's responses have been entered and 


the system is ready, issue the command: 


EDIT ASNF JCL 


Then type in, or read fac, the program which follows. 


// Standard JOB Card (User's Manual Sec. 3.3.2.1) 

// EXEC PGM=!1 EBGENER 

7/SYSPRINT DD SYSOUTSA 

//SYSIN DD DUMMY 

//SYSUT1 OD OSN#F0718.ASDF, VOL=SER=DISKO2,UNIT=3330, 

// DCB=(RECFM2FB, LRECL280, BLKS!I ZE=3200),D1SP2(OLD,KEEP) 
//SYSUT2 DD DSN#F0718. ASDF, VOLSSERSNSXONOL,UNIT#22314, 

// DCB=(RECFM=FB,LRECL280, 8LKSIZE23200),DISP2(NEW, KEEP), 
// SPACE2(TRK, (40,1),RLSE) 

// EXEC PGMslEBGENER 

//SYSPRINT DD SYSOUT2A 

//SYSIN DD DUMMY 

//SYSUT1 DD DSN2F9718.ASDF 191, VOLSSER201SK02, UNI T23330, 
// DCB=(RECFM=F8, LRECL280, 8LKSIZE=800) ,DISP2(OLD, KEEP) 
//SYSUT2 DD DSN2FN718.ASDF191, VOL2SER=0SX001, UNIT#2314, 
// DCB=(RECFM2FB,LRECL#80, BLKS! ZE2-800),91SPaNEW, KEEP), 
// SPACE2(800,(3,1),RLSE) 

// 


Once this program has been correctly read in to the 
user's disk, ft must be submitted to the operating system 
using the standard SUBMIT command: 

SUBMIT ASDF JCL Pl XXXX 
XXXX Is the user's computer center number. If the SUBMIT 
function is working the following messages will be returned: 
we CARNS XFERED TO 0098P ee 
FNAME FTYPE P1 SENT FOR OS BATCH JOB 
After a short delay the system will send a second message: 
FROM 91998P33 : JOB SENT TO OS WITH 21 CARDS WRITTEN (Time) 


The user should then revert back to section B- above which 


Ly deck containing this program ts available from 
Professor Kirk, Sn 302. 


discusses the TOCP command. After giving the system a period 


of time to execute the JCL job which was just submitted, 


indicate that the ASOF 


execution of the TOCP command will 


flles are present. If the operating system has not completed 


executing the transfer, the TOCP request will respond as In 


paragraph two above. 


B. LOADING ASDF 


Once ASDF resides on disk 9SX201, the remaining required 


steps are simple. The user issues the command: 


EXEC ASDF191 


The system will type a series of lines onto the terminal. 


At the conclusion of these typewritten lines the user is 


prompted to issue the next two commands: 


LOGIN 192 P 


The system responds: 


192 REPLACES P (191) 


The next required entry is: 


EXEC ASDF192 


r These two commands will] cause a very long list of CP/CMS 


commands to be typed on the screen. Each of these CP/CMS 


\ commands is automatically executed by the operating system 


in the process of loading the ASDF programs. There are no 


furthrer commands required from the user to complete’ the 


loading process. 


mee 


“hen the entire ASDF program nackage is loaded, the user 
will be oresented with a title page and the Instruction: 
"HIT SPACE AND RETURN TD CONTINUE." Once the user has typed 
"space" and "return" the display wlll write out a listing of 
the commands at his disposal. These commands have already 


been discussed In chapter four. 
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Vi. USING ASDF 


Each of the following sections constitutes an example 
for one of the three input modes. Each example uses oniy 
one input mode. tn sceneral, these modes may be combined 
during the generation of any particular filter. The mode 


selected should be the one most convenient to the user. 


A. EXAMPLE ONE - INTERACTIVE GRAPHICS 

This first example uses the ASDF program in_ the 
interactive sraphics mode. The objective is to observe the 
responses of a digital filter which has two real zeros, one 
at the orlgin and another at +9.8, two complex poles located 
at (9.8 +/- j0.2), and the filter vain is 1.000. The 
interactive graphics mode is used since the exact locations 
of the poles and zeros are not critical. For the purposes of 
this example, It is assumed that the user has already 
followed the login nrocedure described in the ahove 
paragraphs, and has’ loaded the ASOF program. The user has 
already seen the title page, and the two pages of command 
descriptions. The screen is clear except for the top left 
corner which shows: 

#ASDFeRE 


This eroup of letters ts the sfenature for the ASNF proeram 
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and indicates that the program is” ready to execute the 
user's next command. 

Since we must first enter the poles and zeros into the 
z-plane the appropriate command to type is: 

POLZRO 

The ASDF program will then draw the screen as shown in 
figure 6-1, and generate an audible tone. Initfally the 
cursor will be located at some arbitrary nosition or perhaps 
may not appear on the screen at all. Movement of the thumb 
wheels moves the cursor position. Since we wish to execute 
in the interactive graphics mode, to add areal zero, we 
move the cursor (intersection of the two lines tracing 
across the screen) into the box in the left column labeled 
ARZ, (Add a Real Zero) figure 6-2. We then type = any 
keyboard character except "return. The trace cursor will 
disappear, and the flashing alphanumeric cursor will appear. 
The program will print a carat pointing to the box selected 
by the user. The flashing cursor will move into the box 
selected. 

The user then types the command "snace", "e", "return". 
The flashing cursor disappears, and the trace cursor 
reappears, The absence of an audible tone at this time 
indicates that the trace cursor is to be positioned within 


the init circle. 


“The 4012 terminal in the computer center requires a 
"“control-s" after every response. 
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The trace cursor is adjusted to the nosition of the desired 


real zero, (the origin in this case), as shown in figure 


6-3. <Any keyboard character except "return" is then typed. 


The ASDF program plots the zero at the origin, and returns 


cursor with the audible tone. The audible tone 


the trace 


reminds the user that the trace cursor is to be positioned 


Since another real zero is to be 


in the command array. 


entered, the cursor is positioned as in figure 6-1b again, 


the same steps are followed. When the trace cursor 


and 


audible tone, oosition it within the 


reappears without the 


unit circle as near as posstble to the point (0.8 +/- j0.0) 


as shown in figure 6-4. (To enter the zero precisely at 9.8 


The steps for 


the rectangular or polar mode would be used.) 


entering the complex pole pair exactly paralle’ those for 


the previous entries, and are depicted in figures 6-5, and 


6-6. 


the 


zeros have been entered, 


the 


Once all poles and 


trace cursor may be aligned with any command box, and a 


character typed (not “return''), Once the alphanumeric 


cursor enters that box, the user types "space", 


"return". The screen is cleared and the user is asked if he 


wishes to change the default filter gain. Since the default 


value is 1.0, the user responds '"n'"', the ASDF program clears 


and responds with: *ASDF-E 


in the normal is for the user to 


The next step sequence 


look at the ftiter characteristics. To display these 


UNIT CIRCLE COMMANDS 
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Trace Cursor Position Required To Place A Real Zero At The Origin 


Figure 6-3 
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characteristics the command typed is: 
RESPONSE 
To set a hard copy of the filter responses the command: 
EXEC HRDSCPY 
is typed. Figures 6-7a, b, c, d, e, f, and #2 show the 


pertinent portions of the output. 


3. EXAMPLE TWO = RECTANGULA® MODE 

This second filter example is taken from [Ref. 5], page 
222. The filter is a fourth-order low oass Chebyshev. The 
derivation of the transfer function is described in the 
text. The given transfer function is factored to vive the 
destred pole/zero locations. The filter has two complex 
zeros at (-1.0 +/- 39.9) and two complex pole pairs, one 
pair at (.7498 +/- 5.5348), and one pafr at (.7774% +/-~- 
3.2120). The filter tain ts siven as 0.901836. 

In a manner exactly paralleling the previous example, 
the user selects the command ACZRCT (which stands for Add a 
Complex Zero - Rectangular Coordinates), by positioning the 
trace cursor in the appropriate box, as in figure 6-8a, 


Once a keyboard character Is typed, and the alphanumeric 


Ne") and "return" 


cursor moves into the box, the "space", "e 
keys are typed. The ASOF nrogram ‘vill erase the screen and 
draw the display shown in flgure 5-8b. The user then tyoes 
in the rectangular coordinates of either complex zero in 


the first complex zero nalr, figure 6-8. The system will 
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ates 


1 area x TERE ee 


ask whether or not the root location has been correctly 
entered. The user responds 'y' for yes, ‘'n' for no. If the 
root has been incorrectly entered, an 'n' response will 
rewrlte the screen for another attempt. Upon receiving a 
‘y' response ASDF rewrites the unit cirele and command 
array, and plots the new root. The second complex zero is 
entered exactly as the first. 

The complex poles are entered as were the zeros. The 
trace cursor is first placed [fn the box ACPRCT. (Add a 
Complex Pole - Rectangular Coordinates), as in figure 6-9. 
Figure 6-10 shows the unit circle after all poles and zeros 
have been entered. The POLZRO command is terminated by 
aligning the trace cursor with any box, typing a character, 
wafting for the flashing alphanumeric cursor to enter the 
box, and typing "space", "x", The system responds by asking 
about the filter cain. In this example the filter gain is to 
be changed to .001836, so the correct response is 'y'. The 
new filter cain, 9.001836, is typed to terminate the POLZRO 
command. 

To view the filter responses, the user enters the 
command: 

RESPONSE 
when the wASDF-RE signature appears. The hard copy 
responses of the filter are shown in figures 6-lla through 


6-llg. 
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C. EXAMPLE THREE - POLAR MODE 

This example uses the polar mode to tenerate an all cass 
filter. The filter complex pole oairs are to be located at 
(9.2 at +/- 4&5 degrees), and (9.5 at +/- 30 degrees). The 
fllter zeros are to be at (5.0 at +/- 4&5 desrees), and (2.0 
at +/- 30 degrees). The filter zain is to be 0.01. The 
first step is to change the specifications in desrees to 
radians. Thirty degrees equals .5235 radians. Forty five 
degrees is .7854 radians. The entry of moles and zeros 
exactly parallels the above examples, however, the desired 
comands are ACPPLR, (Add a Complex Pole - Polar 
Coordinates), and ACZPLR, (Add a Complex Zero - Polar 
Coordinates). Figure 6-12a shows the trace cursor nosition 
for ACPPL&. Figure 6-12b shows the trace cursor position 
for ACZPLR. The angular vart of any pole or zero entered in 
the polar mode must be in radians. Figure 6-12c shows the 
entry of the first complex zero nair. Figure 6-12d shows 
the entry of the first complex nole pair. Figure 6-12e 
shows the change in filter zaln from 1.00 to 0.91. Figures 
6-l3a through 6-13¢ show the hard copy responses for this 
filter. It should be noted that zeros which olot outside of 
the unit circle are indicated by stars. ATI such zeros are 
plotted at a radius of 1.1 along the same angular component 
as the actual zero location. Zeros are plotted this way for 
two reasons: first to make the output display consistent for 


all relevant filters, and second to emphasize zeros which 
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are close to, but outside, the unit circle. The importance 


of distinguishing between zeros very near the unit circle, 


either interior or exterior, rests on the extreme 
sensitivity of the nhase of the frequency response as a zero 
moves from the interior, onto, and to the exterior of the 


unit circle. 


Vit. ADDITIONAL EXAMPLES 


The purpose of this chanter is threefold. The primary 
objective fs to provide the user with a series of plots 
demonstrating the general form of those tyoes of filters not 
considered in chapter’ six. Additionally, this chapter 
demonstrates the responses sfenerated by poles located at a 
series of positions within the unit circle. It further 
provides examples for user comparison when learning to 
interact with the z-plane. 

The filters considered in the examnples of chapter six 
include a fourth-order low-pass filter and a second-order 
all-pass filter. The sequence of firures 7-la through 7-1le 
shows the general pole zero oplacements for a high-pass 
filter, and associated responses. Similarly, fieure 7-2a 
throush 7-22 show the resnonses for a band-pass filter. 
Fieures 7-3a throush 7-32 show ai band-reject or notch 
filter. 

The secondary purpose of the chapter is to show the user 
the general forns of time and frequency resoonses as the 
filter noles move in the z-plane. The’ seatience of ficures 
7-4 through 7-11 shows the responses for a real nole at 
various locations on the real axis. Fisures 7-12 throuzh 
7-19 similarly show the responses of a complex nole pair 


moving across the z-nlane. Frequency responses are not 
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shown for poles on, or outside of, the init circle since the 
frequency resnonse is not defined for these cases. Figures 
T-17, 7-79, and 7-21 show the chanse in resnonses for a nair 
of complex voles which move awav from the rea? axis. 

The final intent of the chanter is to nrovide the user 
with additional figures for comnmarison should he need 
further practive in manioulating poles and zeros in the 
z-olane. tt is imnortant toa note that the ASS nroeram wilt 
not allow the user to position ovoles on or outside the unit 


circ e. This limitation results in several olots which can 


not be reorodiuced by the user. 
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Figure 7-10b Unit Sample Response For A Single Real Pole At (1.0) 
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VIITt. tNTERPRETATION OF ASDF OUTPUT 


When the user initially confronts the ASDF program [t 
ts expected that he will attempt to duplicate some of the 
examples presented earlier [In this thesis or similar 
examoles found [fn the literature. These tyoes of 
applications of ASDF are straightforward and the results 
zenerated can be compared with published results. 
Interpretation of the output plots is seldom confusing. 

As the user becomes more familiar with the use of the 
ASOF program and the positioning of the poles and zeros in 
the z-plane he will undoubtedly attempt to generate filters 
with sharper transition regions, and greater attenuation In 
the stop band(s). In general these filter Improvements are 
accomplished by increasing the order of the system = and 
posttioning the poles closer to the unit circle. As the 
order of the system increases, and as the poles of the 
filter become closer to the unft circle, the filter becomes 
Increasingly subject to the difficulties associated with 
the finite precision § arltthmetic used in the filter 
Implementation. 

A tacit assumption made very early tn most digital 
signal processing courses (and often quickly forgotten) Is 
that infinite precision arithmetic ts avallable In 


performing the calculations. This assumption breaks down 
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when digital filters are realized in hardware or simulated! 
In general there are three sources of error which may affect 
the performance of digital fliters. Nue to the finite 
precision representation of the coefficients, the values of 
the implemented filter coefficients are only approximately 
equal to the desired values. Second, errors may be generated 
in the arithmetic operations. In general, the product of 
two n-bit numbers results [na 2n-bit result. An error 
occurs when this 2n-bit number’ is converted to n-bits for 
succeeding computations. The third type of error 
encountered Is the analog to digital conversion required for 
processing analog sienals with digital fllters. Only the 
first two of these errors is of concern in the ASDF program. 

White an in depth analysis of finite precision effects 


ts beyond the scope of this work, let us consider the filter 


which has elght zeros located at the origin, and elght poles 


located at -0.9. Entering these locations into the z-plane 
using the ASDF program generates the filter responses [In 
fizures 6-la, b, c. <A superficial glance at these responses 
Indicates that this filter is unstahle. The filter must be 
stable, however, since all of the noles Iie within the unit 
circle. 

There are several phenomena involved [In this example. 
The first is the dynamic ranze limitation of the system on 
which the filter fs implemented. The 1BM-360 ts limftted to 


computing with numbers between 16°" and S160? << Oriee the 
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larger of these limits Is exceeded, the computations are 
terminated and It is not clear whether the designed filter 
Is actually unstable or has merely exceeded the dynamic 
range limitations of the machine. The dedicated user may 
decide to modify the prozrams presented In the appendicies 
to Investigate the time responses of such systems. Clearly 
the form of the time responses will not be significantly 
altered if a modifted unit sample response 

1 WO We My Oe Se: es 
and a modified step response 

1x 19t9 1% 107% 1x 19% 

is used. 

{f these modifications are nerformed, two other 
phenomena become apparent, each relating to the finite 
precision implementation of the digital filter. The first Is 
that, in general, digital filters have deadbands, the second 
that the filter output may degenerate to a limit cycle. The 
original works by Blackman [Ref. 1] . demonstrate that 
digital filters exhibit both of these effects. Extensive 
coveerage of these phenomena is not apvpropriate here, 
however, a dead-band in a dieital filter can be 
simplistically described as a span of flliter Input values 
for which the filter output does not chanre. The limit 
cycle effect is evident when a digital filter responds to a 
steady state input. For constant inputs, the filter output 


may run through a evcle of values. The pronensity for 
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experiencing limit cycles Is obvious in the calculation of 
the step response. A close look at the unit sample response 
shows that after a sufficiently Jlaree number of Input 
values, this unit sample sequence ‘ecomes in reachine 
steady state, a constant. 

When observing the output response of the examole 
filter, there Is still another possibility for explaining 
this output. This explanation revolves around the fact that 
the poles in the Implemented filter are not at the locations 
expected. In fact, the poles may actually move outside the 
unit circle even thouzh they were entered within. In 
general, there are only a_ finite number of pole locations 
oossible with finite word length coefficlent 
representations. The effect of this limitation ifs best 
explained using another example. Consider a fflter with 
three zeros located at 9.5, and 0.2 +/= J0.5, and four poles 
located at -0.999, -0.f, and -0.69 +/- J0.69. The infinite 
precision transfer function for this filter Is expressed in 
equation 8.1. 


2 = .9 22 » 149 2 = .145 (8.1) 


z+ * 2.979 22 # 3.75822 2% + 2.3497398 Z + .57074868 


Now consider the finite precision case where each of the 
coefficients nust be represented by a six bit number. The 


pole at -9.6 Is expressed In binary (exclusive of sin) as 
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.199119N1109. . . and Is truncated to six bits as .100110. 
When all of the coefficients are so represented, the 
transfer function is expressed by equation 8.2. Clearly the 
poles and zeros of this second transfer function are not the 


73 = .890625 2° +.474275 Z = .140625 (8.2) 


24 « 2.9375 2° + 3.75 22 * 2.3125 Z + 5625 


same in equation 8.1. In general the differences will change 
all of the filter responses. In fact for this particular 
case the new poles of the system are at -9.79522 +¢/- 
J0.16677, and -9.76353 +#/- J0.698616; this filter Is 
unstable! 

The reason for the shift In pole locations Is clearly 
portrayed by the filter whose characteristic equation Is 
equation 8.3. 

(z + .995)8 a9 (8,3) 
This filter has eight poles at -0.995. Note what happens to 
the filter when an error ts. Introduced by truncating the 
infinite precision coefficient values to thirty-two bits. 
For simplicity restrict the discussion to the ramifications 
of an error generated in the value of the constant (.995)° . 
This error can be as large as 2x 10 1°, with this error 
explicitly Included, the new characteristic equation will 


be: 


(z + .995)98 © 2% 107t9s (8.4) 


The roots of equation (8.4) are complex, and ite on a circle 
of radius 0.06 around the point -0.995 +/+ J 0.9. Figure 
8-2 shows that three of these roots lie outside the unit 


circle. Thus with a finite preciston implementation of the 


filter (even with 32 bits of onrecitstion), tt fs possible for 


the pole locations to move so far that the fliter actually 
becomes unstable. 

The results of this section as well as four types of 
ASDF outputs which have been observed, are summarized In 
figure 8-3. Trace A shows the unit sample response of a 
stable filter. Trace B shows the response of a filter which 
exceeds the dynamic range limitations of the Implementation 
system. Trace 81 Is for an unstable filter and B2 is for a 
stable filter caught In a limit cycle. These two traces can 
be observed or explored only by modifying the unit sample 
sequence. Trace C shows a stable filter which falls into a 
limit cycle below the dynamic range Itmit of the machine 
upon which the implementation has been performed. Note that 
the performance of a fllter above the dynamic ranze 
limitation of the {BM 369 (or any other machine) cannot be 


observed, unless the ASDF orogram Is modified. 
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1X. SUMMARY AND CONCLUSIONS 


The Advanced Simulation of Digital Filters has been 
implemented on the 18M 360, and Is Intended for’ the use by 
persons beginning their study of digital slenal processing. 
The programs included within the package have been designed 
to make the user's ftnteraction with the z-plane as simple as 
possible. While the programs themselves are available for 
study, all of the computations are invisible to the user, in 
an attempt to develop the user's Intultion of the discrete 
domain without overwhelming him with arithmetic 
computations. 

The ASDF package provides the user with an interactive 
method by which the filter cole and zero locations can be 
manipulated and the effects observed of such movement. ASDF 
can simulate fllters of up to tenth order, and is punctuated 
with numerous comments and Instructions which are to assist 
the user. The [Interface with the z-plane is avallable in a 
rapid interactive eraphics mode, as well as [n_ the slower 
polar and rectangular modes, each providing freater position 
location accuracy. The filter responses generated by ASDF 
are presented asa serles of displays directly on the 
Tektronix 4012 terminal. Versatec plot capability nas been 
integrated into the system to orovide high quality coples of 


the system responses. 
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A. CONCLUSIONS 


There are several points which are worth considering for 
persons Intending to pursue this or other areas of 
interactive systems development with the graphics terminals. 
While comments may apoear obvious they are stated to 
minimize the difficulties which can be encountered. First, a 
serfous attempt should be made to define the entire data 
structure before beginning the programming stage of the 
work. All phases of the effort must be planned before 
commencing with the actual programming to fnsure that the 
data structure selected ts sufficiently general, and 
sufficiently flexible to endure the additions and changes 
whtch Inevitably occur when inittal objectives are 
overlooked or changed. Second, it its Imperative to master 
all of the software tools available and applicable to the 
project. There are features of the PLOT 10 software which 
exfst and were not used In the ASDF program. This 
shortcoming is due to Insufficient a priori knowledge. For 
example, the writing of text oortions of the screen 
presentations should have been much sImpler. Finally, it is 
important to research all of the attributes of the devices 
used for the final problem solution, and direct the solution 


to take advantage of all the system capabilities. 
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B. AREAS FOR FUTURE WORK 


There are several areas in which this orogram could be 
Improved or extended. The possibility exists for 
transferring the ASDF programs to the Tektronix 4082 
computing system. This move would avofd many of the 
shortcommings associated with the current 1BM 360 system. 
The user would be relatively Immune to system crashes, and 
the constraints placed upon him by concerning the 18M 
system's being operational and/or available. There are 
several disadvantages associated with this nove. With the 
advent of the new computer system the majority of the 
current problems concerning system rellabilfty, and 
constraints on available access hours should he eliminated. 
Furthermore, transferring these programs to the 4082 system 
will require significant rewriting to make execution 
possible on a system which has a much smaller virtual 
machine. Furthermore to fully benefit from the 4082 system, 
the interfacing routines should be rewritten to. take 
advantage of the 4082 refresh capability. It seems that all 
of the effort involved for this change is not worth the 
expected improvements. 

The ASDF program is now useful for analyzing digital 
filters. It seems Iltkely that this capability could be 
expanded into a method for computer-aided design of digital 


filters. There are standard transformations which allow the 


fliter desftgner to man analog filter requirements into the 


appropriate digita) requirements, and as well there are 


transformations which make ovossible the computation of 
digital filter transfer functions from known analog transfer 
functions. There are also a variety of optimization 
techniques referenced In the literature which could be 
implemented. 

The other area of obvious expansion for the ASDF program 
is to include the capabiltty for overforming spectral 
analysis. tt would be worthwhile (though perhaps not 
particularly involved) to include in the system the ability 
to calculate Discrete Fourier Transforms, and graphically 
display the results. Additionally, [It would be possible to 
demonstrate the correspondence between performing 
convolution in the time domain directly and performing the 
convolution fn the frequency domain by multiplying the 
DFT's. The relation between the circular and = Iinear 
convolution could also be graphically portrayed. 

As a final precautionary note, it is not possible to 
overestimate the length of time required to senerate well 
documented, correct, useful orograms on an _ “irteractive 
basis. Changes [In system parameters, and the limitatltons in 
system software can make the oprocess of simplifying the 
user's interactive responsibilities a very difficult process 
for anyone not thoroughly acquainted with all avallable 


resources. The myriad of possible fncorrect user responses 


which must be anticipated make the problem extensive. 


APPENDIX A: PROBABLE SOLUTIONS FOR DIFFICULTIES 


SYMPTOM: RESPONSE command returns an error message upon 
execution. 

SOLUTION: Check the virtual machine size with the command 
"CP QF" to tnsure a virtual machine stze of 4OOK. If the 
user has only been allocated 256, logoff and log back on as 
directed tn chapter five. Load ASDF as described in chapter 


six. 


SYMPTOM: User has inadvertently entered the control 
program, te., CP. 
SOLUTION: Type in "bh" and “return," Normal ororram 


execution should resume, 


SYMPTOM: The carat appears with the blinking cursor, but it 
is not located {nn the command array. 
SOLUTION: Move the trace cursor into the command array, and 


type a "space." Normal operation should resume. 


SYMPTOM: No terminal activity. 

SOLUTION: First, check to see whether or not the system has 
crashed, If the system has crashed, logon as [n_ chapter 
five, and load ASDF as in chapter stx. If the system has not 


crashed, move the cursor to the center of the screen and 
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type "space" and "return." tf normal operation has not 
resumed hit the "break" key. When the CP prompt appears, 


type "b.'' Normal operation should resume. 


SYMPTOM: An error Is detected In the execution of ASDF191 
EXEC indicating that no temporary disk space is available. 
SOLUTION: Check with the consultant or a systems 


programmer. 


SYMPTOM: The user can not escape froma delete pole or 
delete zero command in the interactive graphics mode (IAG). 

SOLUTION: Implicit in each request to delete from the 
z-plane is a request, validating the change to he made. In 
the delete IAG mode the trace cursor returns after 
illustrating the root to be deleted. Type "y'' if the correct 
root was tdentified. Type "n'" if the wrong root was 


identified, and attempt the procedure azatn. 


SYMPTOM: The user did not enter the correct amount of data 
and the program terminated on a 217 error. 

SOLUTION: Type "b" and "return." The program should 
restart with the ASDF sfenature. At this point reenter the 


command desired. If executing the POLZ20 command [ft Is 


advisable to execute the list (LST) command to check on the 


current filter status. 


APPENDIX Bs IMPORTANT SUBROUTINES AND VARILABLES 


The subroutines and variables described in this anpendix 


are applicable to all programs within the ASDF vackage. 


A. SUBROUTINE NAMES 
The subroutines within the ASDF package are divided into 
three major catarorfes: 1) Plot-10 software, 2) Versatec 
software, and 3) routines generated explicitely for the ASDF 
programs. 
1. Plot-10 Software Subroutines 
The Plot-10 subroutines are all adaquately described 
in reference 10, and are ltsted below. 
ERASE MOYABS ANMODE 
RECOVR HOME LINEF 
NEWLIN SMENDO VNENDO 
POLTRN MOVEA DRAVSA 
LINTRN DRAM NRVREL 
MOVREL VCURSR BELL 
POINTA SCURSQ INIT 


TABHOR TABYVER SEELYC 


FIN FINITT DRANA 


DRYA8S 


2. Versatec Subroutines 
The Versatec subroutines are listed below, and are 


explained in detail in reference 7. 


PLOTS PLOT WINDOW 
SCALE AXts SYMBOL 
LINE 


3. ASDF Subroutines 

Each of the following subroutines have been 
generated for the ASDF programs. 
CRRMMM - Add and delete poles and zeros from the z-plane. 
The change "C" can be either "A" for add or "D" for delete. 
The roots "RR" can be "RZ" for real zero, "CZ" for comolex 
zero, "RP" for real pole, or "CP" for complex pole. Mode 
"MMM"! can be "IAG" for interactive graphics, "POL" for 
poles, or "RCT" for rectangular. 
FLTR2 - Prints a listing of the pole and zero locations of 
the filter being considered, 
ASUBK - Generates the a, coefficients required to evaluate 


k 
the transfer function. 


BSUBR - Generates the by, coefficients required to evaluate . 


the transfer function. 

ADJUST - !nsures that the correct number of zeros located.at 
the origin are [fncluded in the coefficients ay and b.- 
RSPNSE - Computes the time responses of the filter. 

FRYNCY - Computes the phase and magnitude of the transfer 


function for the prescribed filter. 
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PRPPLT - Scales the ordinate data prior to plotting the 
responses. 

PLTIMP - Plots the unit sample response. 

PLTSTP - Plots the unit step response. 

PLTTRF - Plots the phase of the transfer function as well as 
the magnitude In linear and decibel format. 

UNTCR - Plots the unit circle with the appropriate poles and 
zeros. 

START - Initializes all arrays and required variables. 
FINISH - Checks to Insure that the constructed filter is 
causal, and generates a file to save the pole and zero 
locations. 

ANGLE - Draws the angle sign for polar commands. 

ERROR - Controls oromram flow and error messazes for user 
generated errors. 

LST - Prints a listing of the current pole zero location. 
PLTTBL - Plots the poles and zeros on the unit circle. 
CHCKSZ - Insures that roots are avallable to be deleted and 
that the maximum system order will not be exceeded. 

UPDATE _ Makes changes in the POLZRO array to corresoond to 
user changes In the z-plane. 

LOCATE - Searches the POLZRON table for the root closest to 
that which the user wishes to delete. 

STRRTS = Puts poles and zeros Into the POLZRO array. 


PLOTRT =- Plots a single pole or zero on the unit circle. 
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B. IMPORTANT VARIABLES 

The variables described helow pertain to all of the 
subroutines and programs within the ASDF vackage. 
IROOTS(I) = Keeps track of the number of poles and zeros in 
the system under consideration. IROOTS(1) Is the number of 
real zeros. !ROOTS(2) is the number of complex zero pairs. 
1ROOTS(3) is the number of real poles. JIRON0TS(4) Is the 
number of complex pole pairs. 
POLZRO(!,J) - The locations of the zeros are stored for 
values of | from one to ten. The pole locations are stored 
for values of | from 11 through 20. The real part of the 
rectansular representation of a root location is stored with 
J equal to one. The imaginary part of the rectangular 
represention for a root location is stored with J equal to 
two. The magnitude and anele of the root's polar 
representation are stored for J vatues of three = and four, 
respectively. POLZR0(1,5) holds the root type. 
{ERROR - Keeps track of the generation of an error. Once 
the error has been identified the value is set to one. When 
the error has been exolained to the user the value is reset 
to zero. 
A(t) and BCI) -* Store the values of the a. and b 
coefficients. 
AA(1) and BB(I) = Store the double precision values of the 


ay and b.. coefficients. 


YIC1) + Stores the values of the unit sample response to he 
plotted. 

YS(!) = Stores the values of the unit step response to be 
plotted. 

YP(1) - Stores the values of the phase of the transfer 
function to be plotted. 

YMC(t) = Stores the magnitudes of the transfer function to be 
plotted. 

YMOB(I) - Stores the transfer function magnitudes In 
decibels to be plotted. 

EN(I) = Stores the time values against which the dependent 
variables are plotted. 

1Y!IPTS - Holds the number of unit sample points to be 
plotted. 

{¥YSPTS - Hotds the number of unit step response points to be 
plotted. 

NNEN - ts the order of the denominator plus one. 

NNUM - [s the order of the numerator plus one. 

NORDER - ts the order of the system under consideration. 
This can be elther of the numberator or denominator. 

LAST - Holds the number of pofnts to be plotted. 

FLTRGN ~ Is the value that the user assigns to the filter 
gatn. 

ZEROS(I) = Holds the complex coefficients which will be made 


real to generate the b, coefficients. 
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POLES(!) - Holds the complex coefficients which will be made 
real to generate the a coeffictents. 

The remainder of the variables are either obvious, or 
used merely for control of the program flow, te. do loon 


indicies, etc. 
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APPENDIX C. STRUCTURE FOR ASDF TEXT FILE 


The file ASDF TEXT is a composite file made up of all of 
the required programs and EXEC files for running the ASDF 
program package. These files have been consolidated to 
minimize the time required for transfer from DISKN2 to 


NSX001. The structure of the comnosite file is listed below: 


BREAK)O1 

EXEC file ASOF192 

BREAKNO2 

TEXT file for POLZ20 command 

BREAK103 

TEXT file for RESPONSE command 

BREAKING 

Information printed out by the INT20 command 
BREAKIOS 

VTEC1 JCL - used for EXEC HRNSCPY 

BREAK106 

VTEC2 JCL - used for EXEC HRDSCPY 

BREAK107 

Lables for RESPONSE olots 

BREAK108 

EXEC file to execute the EXEC PHINCH command 


BREAKIN9 
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Intttal file FT91FN01 with no poles or zeros 
3REAKN10 

TEXT fille for generating JO8 cards 
BREAK911 

TEXT file for the INTRO command 
BREAK112 

Information printed by OPTIONS command 
BREAKO13 

Information printed on the thtle nage 
BREAKN14% 

EXEC file to execute HRDSCPY command 
BREAK115 

TEXT file for the OPTIONS command 
BREAKON16 


TEXT file for the software page clear 
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APPENDIX D. JCL FOR EXEC HRNSCPY COMMAND 


The Job Control Language required for the execution of 
the EXEC HRDSCPY command Is provided to eive Insight Into 
the flle manipulation utilized by the ASDF package. Should 
the permanent locations of any of the portions of the ASDF 


package change this section of JCL must be changed also. 


JCL for VTEC1 


J/GO EXEC PGM=THEPLT, REGION=180K 

//STEPLIB DD UNI T#3330, VOL2SER291SK93,DISPSSHR,DSN2S1502. ASDF 
//FT95F001 DO DDNAMESSYSIN 

//FTO6GFNN1 DD SYSOUTSA 

//FTISFO01 0D DONAMEsPLOTPARM 

//PLOTLOG OD SYSOUT#A 

//VECTR2 DOD DSN#&&VECTR2,N1SP2(,PASS),SPACE=(CYL,(5,5)),UNIT#SYSDA 
J/VECTRI DD OSN#&*&VECTR1I,DISP#(, PASS), SPACE=(T&K,(1,1)),UNIT#SYSDA 
//GO.SYSIN DD * 


JCL for VTEC2 


//PLOT EXEC PGHM=!EVMAPP,COND#(4,LT) 
//FT15F901 DD DONAME*PLOTPARM 

//STEPLIB OD DSN#®SYS1.VTECPLOT,DISP#SHR 
//PLOTLOG DD SYSOUT2A 

//VECTR1 DD DISP#(OLD,DELETE) ,DSN#*4A&VECTRI 
//YECTR2 DD DISP#(OLD,DELETE) ,DSN#&%&VECTR2 
//SYSVECTR DD SYSOUT#(A,,5555) 

//NECTTAPE DD DUMMY 

//- 
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